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•  here  holds  for  any  function  of  r  and  r ' ,  and  shows  that,  analogously,  one  correction  • 
term  is  required  per  functional  limit  of  the  triple  Integral  in  each  of  the  deriva¬ 
tives  in  the  curl  expression.  This  generalized  three-dimensional  Leibniz  rule  is 
then  tailored  to  solve  the  specific  problem  of  determining  electromagnetic  fields 
in  the  source  region  by  direct  differentiation  of  the  vector  potential. 

Finally,  in  applying  the  generalized  expression  to  calculate  the  fields,  we 
simply  insert  the  vector  potential  to  evaluate  the  magnetic  field,  and  then  insert 
the  curl  of  the  vector  potential  to  evaluate  the  electric  field.  Since  the  shape 
of  the  principal  volume  can  be  chosen  arbitrarily  we  choose  a  convenient  pillbox/ 
slab  to  simplify  extracting  and  evaluating  the  relevant  correction  terms  in  the 
generalized  Leibniz  expression.  Reducing  the  correction  terms  yields  identically 
the  expressions  for  the  depolarizing  dyadic  previously  calculated  for  an  arbitrary 
as  well  as  for  a  pillbox/slab  volume.  This  process  also  reveals  that  the  Leibniz 
correction  terms  are  directly  related  to  the  depolarizing  dyadic  which  represents 
the  difference  between  Maxwellian  and  cavity-defined  fields. 
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Application  of  a  Generalized  Leibniz  Rule  for 
Calculating  Electromagnetic  Fieids  Within 
Continuous  Source  Regions 


l.  INTRODUCTION 


In  the  following  treatment  of  electromagnetic  fields  in  source  regions  we  calculate  E  and  H 
due  to  a  continuous  current  source  by  direct  differentiation  of  the  vector  potential,  using  a 
principal  volume  approach  similar  to  that  in  Yaghjlan's  article  on  Maxwellian  and  cavity 
fields^  .  This  paper,  however,  presents  an  alternative  differentiation  method  and 
Interpretation  of  the  mathematics  governing  transposition  of  differential  and  Integral 
operators  when  the  Integration  limits  depend  on  the  differentiation  variable.  What  most 
calculus  texts  term  the  "generalized  Leibniz  rule"  pertains  specifically  to  this  subtle  point 
when  differentiating  1-D  integrals.  Therefore,  It  seemed  appropriate  to  employ  a  three- 
dimensional  version  to  rigorously  perform  the  curl  of  a  volume  integral. 

After  searching  the  literature  for  such  a  3-D  version  of  Leibniz'  rule  and  finding  no 
previous  examples  we  derived  the  3-D  generadized  Leibniz  rvile  from  first  principles.  When 
subsequently  applied  to  the  vector  potential  in  a  current  region,  she  3-D  generalized  rule 
produced  the  rigorous  expressions  for  the  "cavity"  as  well  as  "Maxwellian"  electric  and 
magnetic  fields.  One  merit  In  using  this  direct  approach  is  that  it  relates  the  additional  terms 
arising  from  the  3-D  generalized  Leibniz  rule  to  the  difference  between  the  Maxwellian  and 


(Received  for  Publication  1 1  Jan.  1989.) 

^Yaghjlan.  A.D.  (1985)  Maxwellian  and  Cavity  Electromagnetic  Fields  Within  Continuous 
Sovirces,  Am.  J.  Phys.  53:859. 
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cavity  fields,  that  Is.  to  the  source  dyadic  term.  This  study  also  encourages  further 
appUcatloiis  of  the  3-D  Leibniz  rule  derived  here,  wherever  differentiation  is  required  of 
integrals  with  limits  that  depend  on  the  differentiation  variables.  Such  possibilities  lie  in  the 
analagous  electrostatics  problem^  .  and  perhaps  in  the  theory  of  fluid  dynamics. 


2.  DIRECT  DIFFERBN ;  lATION  OF  THE  VECTOR  POTENTIAL  IN  A  CURRENT 
REGION 


FVom  the  egressions  for  the  fields  in  terms  of  the  vector  potential  A  of  a  volume  current 
density  J,  with  the  e'*®*  time  dependence  suppressed,  the  H  field  requires  one  curl  operation, 


H(r)>^^VxA(r)j, 


(1) 


and  the  E  field  two  in  succession^ 


E(r) 

where 


1 


p^Jlr)  -  V  X  V  X  A(r) 


(2) 


Air)  -  u 


u„f 

^0^.0  J 


J(r')  l 


dV- 


V  -  V, 


r  -  r 


(3) 


Though  we  will  only  treat  current  sources  initially,  we  will  generalize  the  results  later  to 
Include  polarization  and  magnetization  souices  as  well.  To  imambiguously  define  the 
"Maxwellian"  fields  inside  the  source  region  we  use  a  limiting  principal  volume,  V5,  in  Eq.  (3) 
to  eliminate  the  singularity  of  the  Green's  function  right  at  the  observation  point.  We  may 
relate  these  mathematically  defined  Maxwellian  fields  to  operationally  defined  "cavity"  fields 
by  cutting  a  hole,  Vc,  in  the  source,  inserting  a  test  charge,  q.  measuring  the  force  exerted  on 
the  test  charge  by  the  now  external  sources,  and  calculating  the  fields,  from 


^  Portis,  Alan  M.  (1978)  Electromagnetic  Fields,  Sources,  and  Media,  John  Wiley  &  Sons.  New 
York,  p.  40-41. 
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E  =  ~q^  •  difference  between  the  Maxwellian  electric  field  and  the  operationally 

defined  cavity  electric  field  Is  contained  In  a  source  dyadic  term,  which  emerges  as  the 
ultimate  product  of  this  analysis,  it  Is  Important  to  realize  that  Is  fixed  during  the 
differentiation  with  respect  to  r  whereas  In  Eq.  (3),  Vg  shrinks  around  the  observation  point, 
moving  with  the  differentiation  variable,  r.  Thus,  V  -  Vg  varies  with  r.  The  mathematical 
repercussion  of  this  Is  apparent  whenever  we  take  the  curl  of  A  (that  is,  differentiate  with 
respect  to  r).  Specifically,  the  limits  of  the  Integral  In  Eq.  (3)  depend  on  the  differentiation 
variable,  which  prevents  a  casual  interchange  of  the  V  and  j  operations. 

To  overcome  the  restriction  on  Interchanging  this  differentiation  and  Integration, 

Yaghjlan  changes  the  Integration  variable  from  r'  to  r"  =  r'  -  r,  to  remove  the  dependence  of  the 
limits  of  Integration  of  the  differentiation  variable,  r,  and  proceeds  to  evaluate  the  fields 
using.  Initially,  a  spherical  principal  volume.^  Since  Leibniz'  rule  already  governs 
differentiation  of  integrals,  an  obvious  alternative  to  Yaghjlan's  change  of  variable  Is  to 
extend  and  apply  the  generalized  Leibniz'  rule  to  the  3-D  differentiation  of  the  vector  potential 
volume  Integral.  This  straightforward  method  retains  both  the  singularity  and  the 
dependence  of  the  limits  of  Integration  on  r,  and  still  )rlelds  Identical  results  for  the  fields. 

References  to  multi-dlmenslonal  forms  of  Leibniz'  rule  are  rare  and  apparently  address 
only  the  case  of  fixed  Integration  limits  rather  than  the  "generalized  Leibniz  rule"  which 
allows  variable  limits  of  integration.  For  example,  Osgood  states,  "In  the  case  of  multiple 
Integrals,  we  assume  that  the  region  of  integration  is  fixed®.  Cases  arise  In  hydromechanics 
In  which  the  region  varies  with  the  parameters,  but  the  treatment  does  tiot  belong  to  the 
elements  of  the  calculus."  A  3-D  rule  is  required  for  the  problem  at  hand,  so  at  this  point  we 
embark  on  a  derivation  of  the  3-D  generalized  Leibniz  rule. 


3.  FORMULATION  OF  3-D  GENERALIZED  LEIBNIZ  RULE 


In  Cartesian  coordinates,  the  curl  operator  can  be  written  as  a  sum  of  partial  derivatives. 


'3K,  8K^- 

XV 

+  V 

+  z 

'  J 

dz  dx 

dx  dy 

K  J 

(4) 


where  K,  In  our  problem.  Is  the  vector  potential  volume  Integral.  We  wish  to  apply  the 
familiar  1-D  Leibniz  rule  to  each  of  the  six  partial  derivatives  above.  Consider  first  the 
generalized  1-D  Leibniz  rule: 


®  Osgood.  (1925)  Advanced  Calculus,  MacMillan,  New  York,  p.  463. 
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sLjW  ajtx) 

ajtx)  aj(x) 


Ekj.  (5)  shows  that  exchanging  the  differential  and  integral  operations  adds  (subtracts)  an 
additional  term  for  each  variable  upper  (lower)  limit.  We  will  find  that  the  curl  of  a  volume 
Integral  will.  In  Its  most  general  form,  require  six  additional  terms  (one  for  each  l*mlt  In  the 
volume  Integral);  these  will  generate  a  total  of  36  additional  terms.  Fortunately,  a  proper 
choice  of  principal  volume  geometiy  will  later  eliminate  most  of  these  terms. 

We  derive  the  3-D  rule  In  similar  fashion  to  the  1-D  rule,  appl)hng  first  the  chain  rule  and 
then  the  fundamental  theorem  of  Integral  calculus.  As  an  example,  we  write  explicitly  one 
term  from  Eki.  (4): 


8z  dz 


F  (r.  r'ldV’, 

X 


where  K  depends  on  the  variables  a^a^  (r.y'.z’),  a^.a^  (r.z'),  a^.a^  (r)  and  r.  The  derivative  of  each 
component  of  K  can  be  expanded,  using  the  chain  rule,  as  a  sum  of  the  derivatives  of  each 
parameter,  holding  the  other  six  constant. 


I 


Ml 


3z  dr  dz 


aj.a^.  ag. 

*4'  *5'  constant 


(6) 


The  fundamental  theorem  of  Integral  calculus  Is  now  applied  to  the  volume  integrals  as 
follows: 


da.  da. 


F^  (  r,  x’,  y',  z'  )  dx'dy'dz'  = 


F^  (  r,  X'.  y'.  aj  )  dy’dz'.  (7) 
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similar  expr**<i‘'<ons  hold  for  32  -  ae  as  well.  If  we  apply  the  chain  rule  to  each  partial 
derivative  'n  Eq.  (4),  and  In  turn  the  fundamental  theorem  to  each  term  in  Eq.  (6),  after 
collecting  terms  we  can  write  the  generalized  3-D  Leibniz  rule  In  a  form  analogous  to  the  1-D 
n:le‘*  ,  with  the  help  of  the  Levt-Clvlta  symbol.  Specifically, 


*1 


+  Vi  ^ 

1  n=l 


+ 1  (-1) 
n»3 


+  I  (-1) 

ns5 


aj 


V  J  wag  ./a, 


r^4 

f 

F^(r.  a^.  x'^.  Xgldx'^dXg 

jJas 

fag 

r 

F^(r.xj.a^.x3)dx,dx^' 

Jaj 

r®2 

I, 

(8) 


Note  that  summation  over  repeated  Indices  Is  implied  and  thus  each  surface  Integral  Is  a 
summation  of  12  additional  terms.  This  Is  the  central  result  which  we  will  use  to  calculate 
the  electric  and  magnetic  fields  In  the  source  region  given  the  vector  potential.  We  emphasize 
that  the  field  calculation  Is  done  to  verify  that  the  square-bracketed  terms  In  Eq.  (8)  give  rise 
to  the  source  dyadic  term. 


4.  TAILORING  THE  3-D  GENERALIZED  LEIBNIZ  RULE 


Returning  to  our  specific  problem  of  evaluating  the  vector  potential  integral,  Eq,  (3),  let  us 
choose  Initially  a  disc-shaped  (pillbox)  prlncl^.al  voluiiie  to  evaluate  (he  surface  Integrals. 
Choosing  a  disc  principal  volume,  rather  than  a  sphere  or  cube,  ii  this  case  Is  advantageous 


*  Kaplan.  (1952),  Advanced  Calculus.  Addlslon  Wesley,  Reading.  MA.  p.  221. 
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because  It  greatly  simplifies  the  3-D  Leibniz  expression  by  eliminating  all  but  2  of  the  possible 
30  additional  terms. 

We  wish  to  evaluate  the  vector  potential  A  by  summing  the  contributions  from  the  surface 
of  the  disc  principal  volume  to  the  edges  of  the  source  region,  as  the  height  of  the  disc.  2e. 

llm 

shrinks  aroimd  the  observation  point  r  (Figure  1).  We  will,  accordingly,  substitute  and  ^  ^ 
for  Vjand  on  the  surface  integrals.  Obviously,  the  outer  limits  of  Integration  are  fixed, 

since  the  source  volume  Is  finite.  Furthermore.  In  Appendix  A  we  prove  that  the  truncated 
disc  can  be  replaced  by  the  thin  slab  tn  Figure  2  where  the  disc  radius  extends  to  the  edges  of 
the  source  region. 


Figure  1.  Pillbox',  or  'Disc'  Principal  Volume. 


Next,  divide  the  volume  of  Integration  Into  the  volumes  above  and  below  the  thin  slab. 
Each  of  these  two  volumes  are  singly  connected  and  thus  we  cam  apply  our  generalized  3-D 
Leibniz  rule  lEkj.  (8)J  to  each  region.  The  limits  In  the  x-  and  y-  directions  are  constant,  and 
the  only  variable  limits  remaining  when  Eq.  (8)  is  appUed  to  each  region  are  the  z' 
integrations.  We  now  integrate  over  the  source  volume  from  z’  =  z■^  e  to  Dtop.  and 
z’  =  z  -  e  to  flbot-  and  from  x',  y'  =  to  flright  where  £2  represents  the  distant  boundary  of  the 
volume  current,  J.  The  only  additional  Leibniz  terms  demanding  consideration  are  those 
corresponding  to  the  z'  limits  of  Integration.  In  particular,  our  bounded  current  region 
(excluding  a  slab  principal  volume)  in  Eq.  (8)  yields: 


Vx 


-  JJJ 


V  X  F  (r.  r')dv 


V  -V. 


v-v. 


(9) 


lim 
e  •  0 


IJ 


Fy  (r.  x'.  y',  z'  =  z  +  £)  dx'dy'  -  |  )  Fy  (r.  x'.  y',  z'  =  z  -  e)  dx'dy' 


+  y 


-  (*•.  X’.  y',  z'  =  z  +  f.)  dx'dy'  +  JJ  (r,  x',  y’.  z'  =  z  -  e)  dx'dy 


where  the  limit  representing  the  surface  of  the  principal  volume  splits  into  two  terms, 

one  for  the  upper  slab  surface  and  one  for  the  lower  slab  surface.  This  expression  Is  particular 

to  our  geometry  but  holds  for  any  sufficiently  well-behaved  vector  function  F. 


5.  CALCULATION  OF  MAXWELLIAN  FIELDS 


To  calculate  the  Maxwellian  magnetic  field  H(r)  from  Eq.  (1)  at  points  within  the  source 
region,  we  simply  differentiate  the  vector  potential  A  In  Eq.  (3)  using  the  3-D  generalized 
Leibniz  rule  (Eq.  (9)1 
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up  to  this  point  we  have  dealt  with  the  problem  of  rigorously  Interchanging  the  curl  and 
Integral  operators  In  Eq.  (10),  but  have  simply  assumed  that  the  curl  and  limit  could  be 
Interchanged.  This  assumption  can  be  proven  valid  by  dividing  the  volume  of  integration 
(V  -  Vj)  into  two  regions,  an  exterior  region  (V  -  Vj)  and  an  annular  region  -  ^5),  where  is 
a  small  but  finite  fixed  volume  enclosing  Vs.  Since  V  -  Vj  does  not  depend  on  8,  the  curl  and 
limit  for  this  part  of  the  volume  Integration  can  be  interchanged  Immediately.  If  Vj  is  made 
small  enough  so  that  J  can  be  expanded  in  a  power  series  over  Vj.  the  annular  volume 
Integration  over  Vj  -  Vjcan  be  performed  and  the  curl  taken  before  or  after  the  limit  Is  taken. 
This  procedure  shows  directly  that  the  result  remains  Independent  of  whether  the  curl  or  Umlt 
Is  taken  first,  thus  proving  that  the  curl  and  limit  can  be  Interchanged  in  Eq.  (10). 

To  evaluate  the  four  surface  Integrals,  we  note  that  If  J  w  (upper)  =  J  w 


8 


Oower)  and  if  J  y 


(upper)  =  J  w 
1^  ■  z  +  e  ” 


(lower),  all  the  surface  Integrals  cancel.  Since  J 


approaches  the  same  value  on  the  top  and  bottom  surfaces  of  the  slab,  as  the  slab  becomes 

,  .  ,  . .  (x  -  x'Y  +  (y  -  y')^  +  (z  -  z'Y 

very  thin,  we  need  only  consider  y  ^  ^  /  \  -i 


V (x  -  x’)^  +  (y  -  y’)^  +  (z  -  z')^ 


Obviously,  for 


z'  =  z  +  e  or  z’  =  z  -  e,  (z  -  z')^  =  e^.  so  vjfj.'  _  ^  +  g  =  ^ .  g.  That  Is,  the  Green's  function  has  the 

same  vadue  on  either  surface  of  the  slab,  and  thus  It  follows  (see  Appendix  A)  that  both  pairs  of 
terms  vanish  in  Eq.  (11)*.  leaving  the  rather  antl-cllmacUc  result  for  the  magnetic  field. 


H(r)  =  J- 


llm 

8-n 


J 


V  X  J  (r’)  V  (|r  -  r'l )  dV 


v-v. 


(12) 


As  Yaghjlan  states,  "The  merit  of  the  rigorous  procedure  does  not  become  obvious  until  It  Is 
applied  a  second  time  to  find  the  Maxwellian  electric  field."  ^  The  electric  field  calculation 
proceeds  from  Eq.  (2)  in  the  same  fashion  as  that  for  the  magnetic  field,  except  we  differentiate 
the  vector  potential  twice:  E(r)  «  V  x  V  x  A(r).  The  first  differentiation,  we  know  fi-om  the 
magnetic  field  calculation,  yields  no  additional  Leibniz  terms,  so  we  may  exchange  the 
differential  and  Integral  signs. 


Vxhm  I  jxudVshm  1  Vx  JiudV. 

8 .0  J  8  .0  J  (13) 

V-Vj  V-Vg 


Dlfierentlatlng  again  we  write. 


Vx 


V  X  fini 

8.0 


J 


JydV 


V  -V. 


=  Vx 


u„  f 

8.0  J 


Vx  JvdV’. 


V  -  V. 


(14) 


where,  as  explained  above,  we  may  Interchange  the  curl  and  limit.  To  evaluate  the  RHS  of 


•  llm  I 

*  for  the  case  where  x  =  x'  and  y  =  y',  J  - _ dS'  becomes  Ik  on  both  surfaces 

S,  ^ 


9 


Eki.  (14),  we  apply  the  same  3-D  Leibniz  rule,  this  time  setting  F  ■  V  x  JV  In  Eq.  (9).  Thus, 


E{r)-^ 


llm 

e-tO 


JJJ 


V-V, 


V  X  V  X  J(r')  v(|r  -  r1 )  dV 


+  X 


JJ  (v  X  Jv)y|^  ^  ^  ^  ^dxdy-  -  IJ  X 


dx’dy' 


iz  «  z  -  e 


(15) 


+  y 


-J*J*  (v  X  dx'dy'  + 

.  S, 

s* 

j/ 


To  evaluate  the  two  pairs  of  Leibniz  surface  Integrals  In  Eq.  (15),  we  first  use  the  following 
simple  vector  identity  to  transform  the  Integrand  F  =  V  x  JV; 


VxJva:Vx\|fJ  =  VyxJ  +  v  VxJ  .  (16) 


Noting  that  V  x  J(r')  =  0,  since  V  operates  only  on  r,  we  can  rewrite  the  surface  integrals  In 
Ek).  (15)  as: 


JJ  (vyx  j)y|^  ^  ^  ^dx'dy'  - 


+  y 


dx'dy'  +  JJ  (V.x 


(17) 


Next,  we  examine  the  x-  and  y-  components  of  the  vector  product  V\|/  x  J: 
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=  (18) 


where  the  components  of  the  gradient  of  the  scalar  Green's  function  can  be  written 


(19) 

X,  =  X 

X2  =  y 

X3  =  Z 


vv=I 


-Iklr  -  r1 


[ik|r  -  r'l  -  1] 


1  X  1 


4rt|r  -  rl 


Xi-x'ij 


Recall  that  In  the  magnetic  field  expression,  the  scalar  Green's  function  V  had  equal 
magnitude  on  the  upper  and  lower  surfaces  of  the  slab  principal  volume.  This  was  due  to  the  z 
dependence  being  quadratic  (that  Is,  (z  -  z')^  =  e2  regardless  of  z').  In  the  case  of  the  electric 
field,  though,  the  factor  of  (z  -  z')  corresponding  to  1  =  3  in  W  is  linear.  The  result,  upon 
substituting  the  components  from  Eq.  (19)  Into  Eq.  (18),  and  subsequently,  Eq.  (18)  Into  Eq. 
(17).  is 


.  lim 
etO 


Y  (x  -  x')  (x'.  y'.  z  +  £)  +  ye  J^(x',  y'.  z  +  €) 


-  y  (x  -  X')  (x',  y'.  z  -  e)  +  ye  J^(x'.  y'.  z  -  e)  ]dS’ 


+  y 


Y  (y  -  y')  (x'.  y'.  z  +  e)  +  ye  Jy(x',  y'.  z  +  e) 


-  Y 


(y  -  y')  Jj,  (x'.  y.  Z  -  e)+  ye  Jy(x'.  y'.  z  - 


(20) 


where  the  shorthand  notation  v  = 


_e^l'-"1[lk  |r-  r’l  -  l] 


has  been  introduced. 


4k  r 


At  first  glance.  It  appears  that  the  pairs  of  terms  containing  y  (x  -  x')  and  y  (y  -  y')  will 
cancel  Immediately  while  the  terms  containing  e  will  add.  Indeed,  the  following  evaluation  of 
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these  Integrals  proves  this  to  be  true.  The  analysis  used  to  evaluate  the  terms  in  Eq.  (20)  is 
similar  to  the  H  field  case,  of  Appendix  B: 

(1)  AMthln  a  continuous  current  distribution,  opposing  cuirents  on  the  upp>er  and  lower 
surfaces  of  a  slab  of  infinitesimal  thickness  c  are  equal.  To  be  certain,  however,  we  expand 
J  (z  +  e)  and  J  (z  -  e)  in  a  Taylor  series  about  z: 


Jj  (z  +  e)  =  Jj(z)  +  e  J’j(z)  +  L-  Jj(z)  +  . .  . 


Jj  (z  -  e)  =  Jj(z)  -  Ej  J'(z)  +  +  •  •  • 


(21) 


Subtracting  the  two  gives  26j'i  (z);  adding  gives  2Ji  (z).  Thus,  in  the  limit  as  £-»0,  all  the  current 
terms  become  zero,  because  they  are  each  multiplied  by  a  factor  e  at  some  point.  Therefore, 
whenever  |W  is  shown  to  be  finite,  we  can  discard  the  product  ejjW. 

(2)  The  gradient  of  the  scalar  Green’s  function  W  yields  slightly  different  integrands  for 
the  X*  and  y-  components  than  for  the  z-  component.  Translating  into  polar  coordinates, 
according  to  Figure  3  we  have: 


2Y(x-x')dS'  =  bm  e 


.271  ,.■»  lkVp’^+  ^ 


poo 

Jo  Jo 


e~'  ■  *’  *  ^  ikp'^cos  0'  dp'dO' 
2n  (p'^  +  e^) 


(22) 


-  Um 
e 


___  j.2.  p 

'’Jo  Jo 


(  .2 

(p  -p] 


and 


'.‘Sj  2T£<iS  =Hge£ 


,lkV  p'^  ♦  e 


ikp'dp'dO' 


271  (p‘^  +  e^) 


•271  «oo  ikVp'^  ♦  , J  . jQ 

liTB  I  1  ®  pdpdG 


Jo  Jo 


27c  [  p’^  +  e*  j 


3/2 
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SIDE  VIEW  OF  SLAB 


0 


TOP  VIEW  OF  SLAB 
SURFACE 


I 


i 


Figure  3.  Polar  Coordinates  for  Evaluating  Surface  Integrals  from  Leibniz  Rule  Application 


Because  we  will  be  taking  the  bmlt  as  e-»0  and  each  of  these  integrals  is  multiplied  by  a  zero 
current  term,  we  can  discard  all  portions  that  are  finite  as  e-*0.  First,  note  the  graphs  of  each 
Integrand  In  Figure  4,  revealing  the  validity  of  the  following  approximations: 

a.  Since  the  singularities  only  occur  as  p  ->o  and  e^O,  we  discard  the  finite  portion  of 
the  curves  above  p'  =  d,  thereby  establishing  d  as  the  upper  limit  on  the  surface 
Integration. 

b.  For  very  small  p'  and  e.  eikVp'^  + 
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Figure  4.  Graphs  Depicting  Infinite  and  Finite  Portions  of  Integrands  Involving  Vy. 


With  these  assumptions,  the  first  of  the  four  Integrals,  depicted  in  Figure  4c,  is  finite  and 
vanishes  in  the  limit  as  e-»0. 

Because  tlie  singularity  in  the  remaining  three  Integrals  is  higher  order,  we  evaluate  these 
Integrals  directly.  First,  we  note  that  the  dB'  Integral  conveniently  yields  a  factor  of  2  7t.  All  p' 
Integrals  are  tabulated,  and  are  evaluated  as  follows: 


In  c 


■  -  hm  e  In  (2d)  -  elne  =  0, 

6*0  e  iC 

fThe  last  term  vanishes  by  L'Hospltal's  rule.) 


.  llm  el  P  ^ 


(a)  e^o 


r 


,  .3/2 


-  llm  e 

€.0 


V? 


+  e 


In  d  +  Vd^  +  e*j  - 
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.d 

1 

d 

r  p'dp’  -  lime 

-  1 

—  .  lim  c 

1  1 

3/2  e.O 

JO  (  ,2  2^ 

P  +e 

_Vp  +6  J 

e^O 

0 

Thus,  only  the  last  Integral  In  Eq.  (22)  survives. 

This  non-zero  result  for  the  fourth  Integral  Is.  In  fact,  the  crucial  discrepancy  arising  from 
the  double  curl  operation  in  Eq.  (2).  It  is  partlcxdar  to  the  slab/disc  principal  volume.  We  will 
quickly  see  that  it  generates  the  correct  source  dyadic  term  L§,  for  a  slab/dlsc.  Inserting  the 
results  of  Eq.  (21)  and  (22)  Into  Eki.  (20),  we  find 


llm 

enO 


{'01  Y  (x  -  X')  -  Ye2jJ  dS’  +y  jj 


y(y-y')  2ej’  -Ye2J 


} 


dS-  ^  (23) 


:x  J^fr)  +y  Jy(r). 


Finally,  we  Insert  Eq.  (23)  Into  the  electric  field  expression. 


E(r)  = 


1(06, 


J(r)  - 

M) 


IJJ 


v-v. 


Vx 


Vy  (|r  -  r’l )  X  J(r’)  dV  -  xJ  (r)  -  yj  (r) 

X  y 


and  combine  terms  to  get  the  final  electric  field  result: 
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Vx 


Vvdr-rDx  J(r‘) 


dV  +  zJ  (r) 
z 


(24) 


where  zJz(r)  is  the  scalar  product  of  the  current.  J*  ,  and  the  pillbox/ slab  source  dyadic. 

I 

L5  >  ZZ  • 

It  Is  desirable,  at  this  point,  to  express  the  E  field  In  terms  of  the  electric  dyadic  Green's 

1 

+  1  Y  (|r  -  r'l ).  To  do  this  we  employ  the  identity 


function.  Gj.(r.  r') 

V  X  X  Jyj  a  J  • 
V%=  -k%.  *■  r’j. 


VV- 


V  and  invoke  the  homogeneous  scalar  wave  equation. 


E(r)  =  ^ 
Icoe, 


.  llm  f 

*  •0  J 


J(r') 


k  I  +  VV 


Vdr-rDdV+zJ  (r) 


v-v. 


(25) 


I  J(r )  • 

1 

_ 1 

i 

L  ^  J 

v-v. 


E(r)  =  IcaHji  J  J(r')  •  Gg.  dV 


V  ( |r  -  r'l)  dV  +  zj^(r)/lti)ep 


zJ^(r) 

lQ)e„ 


v-v. 


Analogously,  we  can  express  H(r)  (Eq.  (12))  In  terms  of  the  magnetic  dyadic  Green's  function. 


^Ya^Jlan,  A.D.  (1980),  Electric  Dyadic  Green's  Fimctlons  In  the  Source  Region,  Proc.  fEEE,  68: 
258. 
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H(r)  = 


Urn 

5^ 


J 

V  - 


(26) 


where 


Gj^  (r.  r')  ■  V  X  \j/l. 


6.  CALCULATION  OF  GENERIC  SOURCE  DYADIC  GIVEN  A  PARTICULAR  SOURCE 
DYADIC 


Having  calculated  the  source  dyadic  corresponding  to  a  designated  principal  volume, 
one  can  derive  a  general  expression  valid  for  an  arbitrary  principal  volume,  Vg.  This  process 
Improves  our  appreciation  of  how  the  Leibniz  rule  surface  Integr  als  relate  to  the  general  form 
of  the  source  dyadic,  Lg,  also  a  surface  Integral.  Restated  mathematically,  we  can  show  that 
our  result  for  the  electric  field,  excluding  a  disc  principal  volume,  is. 


E(r)  a  la# 


llm 
0  6-,0 


J(r')  •  Gg,  (r.  r')  dV  + 


z^r) 

toCo 


(27) 


can  be  rewritten  and  simplified  to  yield  the  general  expression. 


E(r)  .  1CM„ 


j 


V-Vj 


J(r')  •  Gg.  (r,  r')  dV  + 


loe^ 


where 


L 


8 


Sg 


ne 

_ r 


dS 
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and 


r  K-r| 

We  accomplish  this  by  circumscribing  the  disc  principal  volume  with  the  arbitrary  principal 
volume®  (or  vice  versa),  as  shown  in  Figure  5. 


Figure  5.  Arbitrary  Principal  Volume  Circumscribing  Disc-shaped  Principal  Volume. 


sYaghJlan,  A.D.  (1982),  A  Delta-Distribution  Derivation  of  the  Electric  Field  In  the  Source 
Region,  E^trornagnetics,  2:161-167. 
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The  task  Is  to  solve  the  second  Integral. 


Itm 

5 


)dV 


''8 -'^8 


by  considering  the  Integral  over  Just  the  Green’s  function. 


G^dV  = 

6.0  J  6.0  J 


-V3 

dV. 


Again,  we  may  bring  the  J(r')  outside  the  Integral  since  the  region  under  consideration  Is 
small. 

Moving  the  origin  to  r  so  that  |r  -  r'|  =  R  =  jr'l,  and  changing  to  spherical  coordinates,  we  see 
the  unit  dyad  Integral  vanish  as  the  limit  Is  brought  within  the  integral  sign: 


llm 

6.0 


J 

v  - 


.IkR 


4TtR 


dV 


=  lim 

6.0 


I 

V  - 


sin  <t»dRd0di{i 


= iini  r 

6.0  J 


I  ^  sin  (t)d§d6d<t) 

4 


V  -V, 


=  0 


8=  2R 

d5.  =  2 

dR 


To  perform  the  VV  Integral  over  the  primed  volume,  we  use  V’  =  -  V  and  apply  the 

the  outside  and  Inside  surfaces  of  the  shaded  volume  In 


J  V  adV  =  J  G'^adS'  to 

Vj  -  V?  which  contains  no  singularity. 


equation 
Figure  5, 
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_  llm  I  V  f  Vw  I  dv  = 

6^0j  6.0 


{-  j  u’  V\j>dS'  +  r  u  VwdS'  > 

Js,  "  Js”  "  J 


_lkR 


From  Ekj.  (19).  we  know  Vw  s  J —  Ik  -  i  I  e  ,  so  we  rewrite  each  surface  Integral  as  a  sum  and 

4rtR  I  R 


solve  them  all  separately: 


llm  f  u'  VwdS'  =  hm  j  u'  e 
8/)J  n  ^  8.0J  n  r 


e^lk 

4jcR 


dS' 


-  hm  r  G'  e 
6*0 J  n  r 


.IkR 


dS'. 


4kR 


For  the  Integral  over  the  surface  of  Vg,  the  surface  element  more  than  cancels  the  singularity. 


- 11-  f  a*'  V«ds'  =  -  ita  li  f  g's  5^  r"  Mn  <.ded» 

‘  1 J  n  «'04ji  J  " 


,  ,  -  VydS 

6.0  I  n  ^  J 


t 


/  9.<. 


lim_l_  r  u  e 

*  *0  4ji  J  " 


_lkR 


dS 


Taking  the  limit,  the  first  Integral  on  the  RHS  vanishes  with  R,  and  e‘*^  =  1 ,  leaving 

J_  f  ^A’  \  L 

u  -J  dS  ■  Lj 


"  R^ 


.J.  For  the  surface  Integral  over  the  disc  principal  volume,  we  write  dS'  In 


cylindrical  (polar)  coordinates  and  take  the  limit  as  |  .0.  eliminating  the  contribution  from 
the  sides  of  the  cylinder,  as  indicated  in  Figure  6. 
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so  the  p  components  cancel,  leaving  J  e  J  be  over  the  cylinder  ends.  Using  the  ratio 

-1^  -^Jb  ^  . 

- UJ —  s  _L-2 — .  we  deduce  that  b  = - £ - .  Obviously.  u_  *  z.  Also  we  note  that  since 

e^  and  u^^lnt  in  the  same  direction,  the  integrals  over  the  ends  of  either  cylinder  add. 


Figure  7.  Cylindrical  Polar  Coordinate  System  and  Unit  Normal  Directions  Assumed  In 

Deriving  Source  Dyadic,  Lg 


Figure  8.  Side  View  and  Cut  of  Disc  Principal  Volume  Showing  Directions  of  Position 

Unit  Vector  and  Components 


nz — Z 

We  now  have,  after  performing  the  d9  integration,  and  setting  e“"'^  P  +  £ 


=  1,  (these 


Integrals  have  already  been  solved): 


Um 

e-»0 


J  u  °  WdS' =  Ik2z ‘‘*5 

D  n 

S 

S 


J 


p'dp’ 


Um  _ 
ZZ  c  ,0  ^ 


J 


p'dp' 


zz- 


Thus,  returning  to  Eq.  (27),  we  simply  substitute 


for 


Um 

6  .0 


Gg.  dV  which  gives 


E(r)  =  icoPo 


Um 

5.0 


J 


J(r')  •  G  (r,  r')  dV  +  J(r') 


zJ  (r) 

f-  + 

icoe^ 


or 
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B(r)  =  lo)Ho 


J 


v-v. 


(r.  r  )dV  + 


J(r) 


1<BE. 


This  procedure  of  circumscribing  a  particular  principal  volume  with  another  of  arbltraiy 
shape  always  results  in  two  terms:  L5.  and  a  term  to  exactly  cancel  the  source  dyadic 
representing  the  original  principal  volu  ae.  The  complexity  of  our  derivation  of  the  disc 
source  dyadic  is  thus  somewhat  redeemed  by  this  simple  connection  to  the  general  form.  In 
fact,  using  this  procedure,  any  accurate  expression  for  the  electric  field  in  terms  of  the 
particular  source  dyadic  for  a  sphere,  disc,  cube,  ellipsoid,  or  any  other  principal  volume 
shape  can  be  transformed  Into  the  general  formula.  In  terms  of  L5,  which  in  turn  can  be  used 
to  find  the  source  dyadic  for  all  other  principal  volume  shapes®. 


7.  POLARIZATION  AND  MAGNETIZATION  SOURCES 


So  far,  we  have  considered  only  a  source  of  electric  current,  J.  To  include  pwlarlzation  and 
magnetization  sources,  we  note  that.  In  general,  the  "electric  current",  Jg.  Includes  both  current 
and  polarization,  but  the  "magnetic  current",  Jh.  consists  only  of  magnetization,  written 


Jj.  =  J  -  10)P 
=  0  +  IcopM 


We  incorporate  polarization  by  simply  substituting  for  Je;  we  Incorporate  magnetization  by 
summing  the  field  equations  and  their  duals,  (  where  E  *H,  H^E.  and  substituting  for  Jh- 

The  total  field  expressions  for  an  arbitrary  principal  volume,  taking  Into  account  current, 
magnetization,  and  polarization,  are  thus; 
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E(p)  =  f  k  .  .  M  .  G^]  dV  . 

J  0 

V-Vj 

»'•  -  IZ  f  [°H  •  ■’e  *  •  m]  dV  -  4  .  M. 

V-Vg 


Having  derived  the  fields  using  the  3-D  Leibniz  rule  It  Is  now  easy  to  appreciate  how  the 
Leibniz  additional  terms  form  the  source  dyadic.  L5.  Comparing  Eq.  (23)  with  Eq.  (25), 


=  J(r)  -  Hm  J  j j  Y  e  (X-.  y-.  z)  dS*  -  j?  J j  y  e  (x'.  y.  z)  dS'.  ^^g) 


On  the  RHS  are  two  additional  Leibniz  surface  Integrals  which  comjjensate  for  the  r 
dependence  of  V5. 

If  we  had  chosen  In  the  beginning  to  evaluate  E  and  H  using  a  cavity  volume,  .  to  exclude 

Ls  •  J 

the  observation  point,  the  source  dyadic  would  never  result.  Hence,  the  source  dyadic  —  is 

Ls  *  J 

the  difference  between  the  Maxwellian  and  cavity  fields,  that  Is,  Em  -  Ec  = - •  It  Is  easy  to 

icoco 

see  why  this  occurs  by  considering  how  the  equation  above  would  read  for  a  cavity  volume 
rather  than  a  principal  volume.  Since  the  cavity  Is  formed  by  physically  removing  a  small 
volume  of  ciurent  around  the  observation  point  r,  J(r)  =  0.  Because  the  location  of  the  cavity 
Is  fixed  during  the  differentiation  of  A,  Vc  does  not  depend  on  r.  Obviously  then,  the 
Integration  limit,  V  -  Vc  does  not  depend  on  the  differentiation  variable,  r.  Thus,  the  regular 
Leibniz  rule  permits  Interchanging  the  differentiation  and  Integration.  Since  the  "generalized 
Leibniz  rule”,  which  Involves  variable  Integration  limits.  Is  not  necessauy  for  evaluating 

Ls  •  J 


cavity  fields,  no  additional  terms  [the  surface  integrals  in  Eq.(28)]  arise,  and 


Icoeo 


=  J(r)  =  0. 


8.  CONCLUDING  REMARKS 


Using  a  principal  volume  approach,  we  have  made  a  thorough  study  of  how 
electromagnetic  fields  behave  In  the  source  region,  by  deriving  and  appljdng  a  generalized  3-D 
Leibniz  rule  to  correctly  differentiate  the  vector  potential.  Section  1  began  by  considering 
only  electric  current  sources  and  emphasized  that  the  difference  between  cavity  and 
Maxwellian  fields  hinges  on  rules  restricting  Interchanging  differential  and  Integial 


25 


operations.  Then,  a  3-D  version  of  the  generalized  Leibniz  rule  was  derived,  to  properly 
perform  the  Interchange  and  was  subsequently  employed  along  with  a  disc  principal  volume  to 
calculate  both  the  magnetic  and  electric  fields.  In  Section  5,  the  fields  for  an  arbitrary 
principal  volume  were  derived  and  found  to  be  identical  to  those  derived  previously,  using 
other  techniques.  Finally,  the  duality  of  Maxwell's  equations  was  Invoked  to  obtain  the  field 
expressions  which  Include  pxjlarlzation  and  magnetization  sources.  In  addition  to  current. 

While  this  exercise  reproduces  previously  known  results  for  the  fields.  In  doing  so,  it  not 
only  provides  a  clear  and  completely  rigorous,  straight-forward  treatment  of  the  problem,  but 
an  intuitive  eiqilanatlon  as  well,  from  both  a  mathematical  and  a  pl^slcal  viewpoint. 
Primarily,  the  extension  to  three  dimensions  demonstrates  the  breadth  of  Leibniz' 
fundamental  rule,  providing  a  glimpse  of  Its  potential  to  enhance  precision  in  volume  Integral 
formulations. 
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Appendix  A 


Proof  of  Validity  of  Exchanging  the  Pillbox  Principal  Volume 

for  a  Slab  Principal  Volume 


Working  In  rectangular  coordinates,  the  most  logical  principal  volume  shape  for  excluding 
a  singularity  of  the  scalar  Green's  funcUoii  y  Is  the  disc,  or  pillbox,  shown  In  Figure  A1  below. 
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To  specify  wiilch  of  the  36  additional  terms  in  the  generalized  3-D  Leibniz  rule  actually 
contribute  to  the  fields,  it  is  advantageous  to  extend  the  x-  and  y-  dimensions  of  the  disc, 
forming  a  slab  principal  volume.  We  wish  to  show  that  the  disc  and  slab  can  be  used 
Interchangeably  as  the  limiting  Vg  in  the  vector  potential  Integral.  To  do  this  we  simply  show 
that  the  vector  potential,  integrated  over  the  shaded  region,  vanishes  as  the  thickness,  2e, 
approaches  zero: 


k  (r) 

shaded 


llm 

5hO 


r'l )  J  (*")  dV  =  0 


(Al) 


or.  In  cylindrical,  coordinates. 


A  (r) 

shaded 


17^ 


J  (r‘)  p'dp’dB'dz' 


(A2) 


where  we  have  moved  the  origin  to  r.  By  definition  of  the  disc  geometry,  .  0  so  we  have 

5 

replaced  5  >0  with  e  «0.  We  assume  J  is  approximately  constant  over  a  thin  Az,  and  perform  the 
z  Integral  first,  noting  that  the  Integrand  is  an  even  function,  symmetric  about  z' s  0: 


llm 
e  >0 


d0’J(r) 


^IkVr'^  +  r-^ 

/  ,2"  ,2  ' 
V  p  +z 


dz' 


Using  the  Inequality 


/•*  IkVr’^  +  z'^ 

f  ®  dz- 

5  f" 

1  «  r.2  .2 

Jz  =  0  Vp  +z 

1  «  /  -2  .2 
Jz  *  0  Vp  +z 

we  can  evaluate  the  RHS  as: 


(A3) 


(A4) 
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Z'  =  O 


(A5) 


llm 

e->0 


J(r') 


')  Infz' 

V 


+  'V  p  +  z 


=  Um  J(r') 


In 


=  J(r’)  [In  (p  )  -  In  (p  )]  =  0. 


Therefore,  Ashaded  =  0  proving  that  It  Is,  In  fact,  valid  to  use  the  disc  and  slab  principal 
volumes  Interchangeably. 


I 
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APPENDIX  B 

Proof  That  ijT  Terms  Cancel  if  They  Are  Equal  and  Opposite, 

Despite  Singularity  in  \y 


Recall  that  In  the  Maxwellian  H  field  calculation,  we  arrive  at  the  following  expression 
after  considering  the  3-D  Leibniz  rule; 


H(r)  =  i 


V  X  J(r')\|f(|r-r’|)  d^r' 


(Bl) 


Where  and  are  the  values  of  V|r|  I  r-r'  I )  on  the  upper  and  lower  surfaces  of  the  slab 
principal  volume. 

The  functions  Vtop  Vb^t  are  equal,  (Vt^p  =  \|tbot  =  y).  yet  still  singular: 
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'»'top  =  v! 


^  |  y(z^^Z■^e) 

V  (x-x)^  +  (y-y)^  +  l-e)^ 


y(z’=z-e) 


'^bot”V 


glkV(x-x')=^+  ty-y  )2  > 
K“*'®  VCx-xf  +  (y-yf  +e^ 


(B2) 


We  wish  to  prove  that  In  spite  of  this  singularity,  the  quantity  In  square  brackets  In  equation 
(Bl)  vanishes.  Since  =  Vbot.  we  write 


llm 

e-*0 


'^ybot 


-  J 


ytop 


dx'dy'  +  y 


i 


J-4  -  J  W... 

xtop  xbot 


dx'dy' 


(B3) 


We  expand  (z'=z  +  e)  and  Jy  (z'=z  ±  e)  in  a  Taylor  series: 


2 

J.  (z'=z+e)  =  J(z)  +  eJ'(z)  +  —  J"(z)  +  ... 

2! 

(z'»az-e)  =  J(z)  -  eJ'(z)  +  —  J"(z)  +  ... 

2f 


Subtracting,  we  have 

J  .  .  -  J  ^  =  -2e  J  '(z) 

ybot  ytop  z 

-J  .  .  =  2eJ '(z). 

xtop  xbot  z 


Substituting  these  expressions  Into  (B3),  we  obtain  for  the  bracketed  quantity 


Urn 

e-K) 


2eJ  ’(z) 

Z 


ydx'dy'  +  y  IJ  ydx'dy'] 
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To  Integrate  over  y,  we  express  the  integrand  in  polar  coordinates  and  Impose  an  upper  limit  d 
on  p’,  since  y  is  singular  for  only  small  values  of  p'  (thus;  we  may  approximate 
as  1). 


JJ 


ydx'dy' 


V  p  +  € 


p'dp'd©' 


-P^  =  2xVd'^ 

~72-  7 

p  +  £ 


(B4) 


Now.  we  take  the  limit  to  show  that  the  bracketed  quantity  is  zero: 


llm 
e— >0 


✓s 

y 


27tV  p'^  + 


~7. 

e 


=  0 


Thus,  Eq.  (Bl)  becomes  simply. 


H(rl  = 


X 

Ho 


llm 


1. 


I  V  X  J(r')y{lr-r'l)dV'.  ■ 
V-V5 
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MISSION 


Rome  Air  Development  Center 


RADC  plans  and  executes  research,  development,  test  and 
selected  acquisition  programs  in  support  of  Command,  Control, 
Communications  and  Intelligence  (C*I)  activities.  Technical  and 
engineering  support  mthin  areas  of  competence  is  provided  to 
BSD  Program  Offices  (POs)  and  other  BSD  elements  to 
perform  effective  acquisition  of  CPI  systems.  The  areas  of 
technical  competence  include  communications,  command  and 
control,  battle  management  information  processing,  surveillance 
sensors,  intelligence  data  collection  and  handling,  solid  state 
sciences,  electromagnetics,  and  propagation,  and  electronic 
reliability /maintainability  and  compatibility. 


